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A polymer placed in chaotic flow with large mean shear tumbles, making a-periodic flips. 
We describe the statistics of angular orientation, as well as of tumbling time (separating 
two subsequent flips) of polymers in this flow. The probability distribution function 
(PDF) of the polymer orientation is peaked around a shear-preferred direction. The tails 
of this angular PDF are algebraic. The PDF of the tumbling time, r, has a maximum 
at the value estimated as inverse Lyapunov exponent of the flow. This PDF shows an 
exponential tail for large r and a small-T tail determined by the simultaneous statistics 
of velocity PDF. 



Introduction. With the development of novel optical methods a number of high qual- 
ity experimental observations resolving individual polymer (e.g. DNA molecule) dynam- 
ics for elongational and shear regular flows have been reported by Perkins et al. (1997)[ 



Smith fc Chu (1998)|and by|Smith et al. (1999)|and|Hur et al. (2001)| respectively. Ex- 



perimental explorations by|Smith et al. (1999)]and the subsequent theoretical/numerical 



study by Hur et al. (2000)| of the shear flow setting had focused on the analysis of the 
power spectral density and the simultaneous statistics of polymer extension, for the 
case in which fluctuations are driven by thermal noise. Additionally, one noticeable 



phenomenon, called tumbling, was reported by Smith et al. (1999) for the regime of 
strong shear flow: The molecule which spends most of the time being oriented along 
the direction dictated by the shear sometimes and suddenly (a-periodically) flips. In 
another experimental breakthrough, a chaotic flow state called by the authors "elas- 
tic turbulence" was observed for dilute polymer solutions by Groisman & Steinberg 
(2000,2001,2004). This flow consists of regular (shear-like) and chaotic components, the 
latter being weaker. Resolving an individual polymer in this steady chaotic flow was more 
challenging than in the regular flow experiment but still it appeared to be an accessible 
task for Gerashchenko et al. (2004) The coil-stretch transition, predicted by Lumley 



(1969,1973) (see also Balkovsky et al. (2000,2001) and |Chertkov (2000)|), was observed 
for the first time, in direct single-polymer measurements by Gerashchenko et al. (2004) 
The statistics of the polymer extension and of the tumbling time were also tested in the 
elastic turbulence experiments, however the phenomenon has not yet been fully explored. 

In this letter, we discuss the statistics of polymers placed in a chaotic flow with a 
relatively large mean shear, that is the flow of the type correspondent to the elastic 
turbulence experiments by Groisman & Steinberg (2000,2001,2004). We assume that 
the effect of velocity fluctuations is stronger than that related to thermal noise and that 
polymers are essentially elongated due to the fluctuations so that the polymer orientation 
is well defined. The main body of the orientational fluctuations occur in a neighborhood 
of a special direction preferred by the shear. Sometimes these typical fluctuations around 
the preferred direction are interrupted by flips, in which the polymer orientation is re- 
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versed. The task of this study is to describe the statistics of the angular orientation and 
tumbling time. 

The structure of this letter is as follows. We begin by introducing the basic dumb- 
bell-like equation governing the dynamics of the polymer end-to-end vector in a non- 
homogeneous flow. If the effect of thermal fluctuations is negligible, the angular part of 
the polymer dynamics decouples from its extensional counterpart and can be examined 
separately. It is convenient to count the angular degrees of freedom <f> and 9 (for in-plane 
and off-plane orientations, respectively) from the direction prescribed by the shear. We 
show that the angular PDF is peaked at some small angle </>, estimated by the average 
value 4>t = ((f)), 4>t <C 1. The widths (in both angles) of the main part of the PDF are 
of the order of 4>t- Then we demonstrate that tails of the joint PDF are algebraic at 
<j>t -C (f), 6 <C 1. We find that this algebraic tail of the individual PDF of <fi is related 
to purely deterministic (i.e. shear driven) dynamics: V((j>) oc sin -2 0. The tail of the 
9 PDF has two competing contributions, one related to deterministic dynamics, oc 9~ 2 
for <§C |6*| 1, and the other one related to stochastic dynamics, oc 9~ a , where 
a is a number, dependent on details of the velocity fluctuation statistics. Then we 
examine the statistics of the tumbling time, r, that is defined as the time between two 
subsequent flips of the polymer. The PDF of t is peaked at a time estimated by the 
inverse Lyapunov exponent of the flow, r t — A -1 . The long time, t > t ( , tail of the 
PDF is exponential, ln[P(r)] ~ — r/r t . The statistics of small tumbling times is related 
to the simultaneous PDF of the velocity gradients. To derive these results we explore 
the close relation between the stochastic dynamics of 9 and <fi and Lagrangian dynamics 
of the flow. separation in the flow. We conclude by discussing the applicability conditions 
for our approach and the validity of the assumptions made in this letter. 

Model. We consider a single polymer molecule which is advected by a chaotic/ turbulent 
flow. We assume that the velocity correlation length is much larger than the size of the 
polymer. (Note, that this condition is always satisfied in elastic turbulence simply be- 
cause the velocity correlation length coincides with the overall size of the flow/apparatus 
that is the biggest scale in the problem.) Then the polymer can be viewed as a material 
point moving along a Lagrangian trajectory. In addition, the polymer is stretched due 
to the flow inhomogeneity. The polymer stretching can be characterized by its end-to- 
end separation vector, R. The stochastic dynamics of the vector R can be examined 



in the framework of the following dumb-bell-like equation (see e.g. Hinch (1977) and 
Bird et al. (1987) ) 



d t Ri = RjVjVi-iRi + Ci, (0.1) 

where 7 is the polymer relaxation rate (dependent on _R, the absolute value of the vector 
-R), the velocity gradient VjVi is taken at the molecule position, and Q is the Langevin 
force. The velocity difference between the end points of the polymer is approximated 
in Eq. (|0.1|l by the first term of its Taylor expansion in the end-to-end vector. This is 
justified by the smallness of the polymer extension R in comparison with the velocity 
correlation length. 

We focus on the situation for which the effect of velocity fluctuations is stronger than 
the effect of thermal noise. Then the Langevin force in Eq. (|0.1|) can be neglected. 
In this case the polymer angular (orientational) dynamics described by the unit vector 
n = R/R decouples from the dynamics of the end-to-end polymer length R and one 
derives from Eq. I|0.1(l a closed equation for n: 



dtrit = rij(5ii - ritn^VjVi 



(0.2) 
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Figure 1. Schematic figure explaining polymer orientation geometry. 

Note that Eq. Il().2[l coincides with the dynamics of a rod-like micro-object immersed in 
the same flow. 

Let us choose the coordinate frame associated with the mean shear velocity, as shown 
in Fig. ^ In this reference frame the mean flow is characterized by the velocity (sy, 0, 0) 
where s is assumed to be positive. The polymer orientation is conveniently parameterized 
in terms of the angles <f> and 9: n x — cos 9 cos 4>, n y = cos#sin</>, n z = sin#. Then, Eq. 
I|U.2I) becomes 

d t (j> = -s sin 2 <f> + ^ , (0.3) 
dtO = — s sin <f> cos <f> sin 9 cos 9 + £g , (0-4) 

where £^ and £g are random variables related to the fluctuation component of the velocity 
gradient. 

Angular statistics. The statistics of the velocity fluctuations is assumed to be ho- 
mogeneous in time. In a statistically stationary velocity field, the angular statistics is 
stationary as well, being characterized by the joint PDF, V{4>,&), which is a periodic 
function of the angles with the period 7r for both <f> and 9. Thus, it is sufficient to con- 
sider V(<p,9) within the following bounded domain (torus), — tt/2 <</>,#< tt/2. We 
normalize the PDF using J d(f>d9V((j>,9) = 1, where the integral is taken over the do- 
main. Note, that according to the structure of Eqs. (|0.3IU.4ll . V(<f>, 9) is symmetric with 
respect to 9 but it is not symmetric with respect to cf>. Therefore, the average value of <j>, 
4>t = (0); is non-zero. In our setting, 4>t is positive. The value of 4>t can be estimated by 
balancing the deterministic and stochastic terms on the right hand side of Eq. I|(J.3|) . The 
weakness of the random term in comparison with s implies <fi t 1. The same quantity 
4>t estimates typical fluctuations of <ft about its mean value. Once one assumes that the 
random terms in Eqs. (|0.3I0.4|> are comparable, it immediately follows that typical value 
of 9 fluctuations is estimated by 4>t as well. 

Note that the equation, d t Ti = rjVjV i: describing dynamics of separation between two 
neighboring fluid particles (moving along nearby Langangian trajectories), leads to the 
same dynamics for the unit vector r jr as determined by Eq. I|0.2|) and, consequently, to 
the same angular dynamics as described by Eqs. I|0.3I0.4[) . For r (the absolute value of r) 
one derives: dt lnr = scos 2 9 cos <j) snv 4> + £|b where £|| represents the direct (as opposed 
to indirect through fluctuations in 0) effect of velocity fluctuations. It follows from Eq. 
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1)0. 3|1 that for typical fluctuations (when <f>, 6 <C 1) competes with sip 2 . Assuming 
that £|| ~ £0 one finds that £|| is negligible in comparison with s<p, because the relevant 
values of 4> are small, 4> <C 1. Therefore, for small (j>,0, one arrives at 9jlnr = s<p. This 
equation establishes the relation between the angular dynamics of the polymer and the 
dynamics of Lagrangian separation. For the Lyapunov exponent, A, defined as the mean 
logarithmic rate of divergence of Lagrangian trajectories, one finds A = {d t lnr) = s<f> t . 

It is natural to expect that the Lagrangian velocity correlation time is r t = A -1 = 
{s4>t)~ l , that is also characteristic time of the and fluctuations. Then, comparing 
the left hand sides of Eqs. 1)0.310.4)1 with the first terms on their right hand sides (for 
<f>, 9 <C 1), one concludes that the angular correlation time can be estimated by the same 
quantity r t . Next, equating the terms on the right hand sides of Eqs. HO.3IO.4fl . one 
derives £0 ~ £g ~ s<j> 2 <C s. The last inequality reflects the assumed weakness of the 
velocity gradient fluctuations compared to the shear rate, s. 

Tails of the angular PDFs. Let us consider the domain \<p\, \9\ 3> <pt, where the 
random terms in Eqs. 1)0.310.4)1 . £0 and £0, are negligible. The angular dynamics is 
purely deterministic in this domain leading to the following dependence of the angles on 
time t 

cot(j) = s(t — to) , tan 9 = c ■ sin <f> , (0-5) 

where to and c are some constants. According to Eq. 1)0.5)1 . the vector n reverses its 
direction as t increases. Therefore, Eq. (|0.5|) describes a single flip of the polymer. 
Due to the assumed homogeneity in time of the velocity statistics, to is homogeneously 
distributed. Recalculating the measure dto into the PDF of the angles in accordance 
with Eq. I)0.5)l . one derives 

v^e) = m 7l'^ ■ (0-6) 

sin <p cos z 

The function U reflects possible variations in c (its statistics), which should be deter- 
mined from the initial conditions for deterministic evolution. These conditions have to 
be found from matching Eq. (|U.5|I to those defined for the stochastic domain \(j>\, \6\ < <j> t . 
One concludes that the function U is sensitive to the angular dynamics in the stochastic 
domain and, respectively, to details of velocity fluctuations, i.e. the function is nonuni- 



versal. Note, that Eq. 1)1). 6fl is identical to the one found in |Hinch fc Leal (1972) in the 
context of a solid rod tumbling in shear flow caused by thermal (Langevin) fluctuations. 

Eq. (|0.5|) shows that in the deterministic regime the angle (f> decreases uniformly 
with time (that is except for the jump from — tt/2 to +tt/2 at t — to). Therefore, the 
stationary PDF for the angular degrees of freedom, V(<t),0), corresponds to a non-zero 
probability flux from positive to negative <f> related to a preferred (clock-wise) direction 
of the polymer's rotations in the X — Y plane. Formally, the probability flux goes out 
through <p = and the same flux comes back (enters) through <f> = tt/2 (tt/2 and 

—tt/2 are identical by our construction) thus keeping the total probability equal to unity. 

The PDF of cf>, P$, can be obtained from the joint PDF: P, p = / d6V((f>, 9). Integrating 
the right-hand side of Eq. (|0.6fl over 9 one obtains the following expression for the tail, 
valid for |</>| 3> <f> t : 

P tjj = J d9 V{<t>, 9) = C(j) t sin' 2 , (0.7) 

where the constant C is of order unity. Let us reiterate that, thinking dynamically, Eq. 
(|0.7(l originates from the deterministic flips bringing from its most probable domain 
~ (fit to the observation angle. Eq. 1)0.7(1 describes the aforementioned probability flux: 
as determined by Eq. 1)0. 3J1 . P^dt<p is constant in the deterministic region. 
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Consider the PDF of 9, P e = J d4>V{<f>, 9). The naive result for the PDF tail following 
from Eq. (|0.6|l is Pg oc \9\~ 2 , provided 1 > fl > (4. However, one should be careful, since 
the expression does not cover a special angular domain, characterized by \<p\ < (f> t 
and \9\ 3> <j>t, which should be analyzed separately. In this domain, one can neglect 
in Eq. I|0.4fl . Assuming also \9\ -C 1, one arrives at <9f ln((9) = —s<f>. In this case, s<f> can 
be treated as a random variable independent of 9 and the above equation leads to an 
algebraic tail, V oc |#| -a , where the exponent a is a positive number of order unity. The 
value of a is sensitive to the statistics of the 4> fluctuations. (Therefore, a is not universal.) 
Let us explain the origin of the algebraic dependence. The algebraic contribution to the 
angular PDF is related to the long (compared to the correlation time Tt) period when (j> 
fluctuates around some negative value, ~ —4>t- (These fluctuations can be interrupted by 
flips.) Then 9 at the end of the T-long period is estimated according to \n(9/<fr t ) ~ Tsfa- 
The probability W to observe such a long a-typical period is estimated by In W ~ — T/r t . 
These estimates, recalculated in the PDF of 9, Pe — dW/d9, give the aforementioned 
algebraic tail. Note that this algebraic tail which originates from the long-time dynamics 
is analogous to the algebraic tail of the polymer extension PDF discussed by Balkovsky 
et al. (2000,2001) and |Chertkov (2000)| 

Therefore, one finds that there exist two different contributions to the PDF tail: one 
related to the deterministic motion, described by Eq. I|H.5|I . while the other is associated 
with the stochastic evolution in the domain, \<fi\ < <j>t,\0\ 3> <j>t- For 1 3> \9\ 3> (f>t, 
both contributions are algebraic, oc \9\~ 2 and oc \9\~ a , respectively. The deterministic 
contribution, oc \9\~ 2 , dominates if a > 2, while the stochastic contribution, oc \9\~ a , 
dominates otherwise. 

Tumbling time statistics. As seen from the expression l|U.5fl . the deterministic pro- 
cess, which actually defines the polymer turn (because <f> changes essentially only during 
deterministic part of the dynamics), is faster than the stochastic wandering taking place 
at small angles, \<j>\, \9\ < (fit- Therefore it is convenient to define the tumbling time, r, 
as the time separating two subsequent crossings in 4> of the special angle ±7r/2, in the 
middle of the deterministic domain. Since the major contribution to r originates from 
the stochastic wandering in the </> t -narrow vicinity of <fi = 0, the position of the r-PDF 
maximum and its width are both estimated by the correlation time Tt = (s</>t) _1 , because 
this is the only relevant characteristic time of the stochastic angular evolution. 

Being interested in the PDF tail, for r 3> r t , one observes that if a flip does not occur 
for a long time, then this delay can be interpreted in terms of the large number, T/r t , 
of independent unsuccessful attempts to pass (clock-wise in <J)) the stochastic domain 
\(f>\ < 4>t- The probability of the delayed flip is given by the product of the probabilities 
of these t/tj events, resulting in the exponential tail of the PDF of r for r 3> r t , 
lnP T ~-T/T { . 

The left, r n, tail of the tumbling time PDF is non- universal because it is sensitive to 
details of the velocity field statistics. Indeed, it is determined by the special configurations 
of the velocity field that force <j> to drift through the stochastic region a- typically fast. 
(Those configurations are vortices with clock- wise rotation of the fluid in the X — Y plane 
leading to negative values of £j, that are larger than sepi in absolute value.) Analyzing the 
anomalously fast revolutions of the polymer, one finds that £^ from the right hand side of 
Eq. I|0.3[) may be considered as time-independent. (Here again, the natural assumption 
is that the correlation time of the velocity field fluctuations is of the same order as the 
inverse Laypunov exponent in the flow, i.e. ~ r t .) Then the direct solution of Eq. I|0.3fl 
gives r — n/ y|^|s, where we assumed that the major contribution in r comes from the 
domain of small (p, (j> <C 1. This estimate holds if s 3> |£<jj| 3> s<$. For 1/s <C r <C n one 
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arrives at the following expression for the PDF of r: 

Pr = ^ (-4-) • (0-B) 

where is the single-time PDF of £4,. 

Conclusions. Let us discuss applicability conditions of our results. One of our as- 
sumptions was that the mean flow can be approximated by a perfect shear flow, whereas 
in reality flow parameters vary along the Lagrangian trajectory demonstrating essential 
spatial inhomogeneities. Even though these variations were not included in our deriva- 
tions, our results remain valid if the variations along the Lagrangian trajectory occur 
on time scales larger than r t and also if the local flow does not deviate strongly from 
the shear configuration. Then, the PDFs discussed in this paper adjust adiabatically to 
the current values of the parameters and become, consequently, slow functions of spatial 
position. If the regular part of the flow is elongational, polymer flips become forbidden 
in the ideally deterministic regime while fluctuations will still generate some tumbling. 
Analysis of the tumbling in elongational, and, generically in any other random flow, will 
be discussed elsewhere. 

We have focused on discussing the dynamics and the statistics of the polymer's angular 
degrees of freedom. However, our theoretical scheme can be naturally extended to also 
include polymer extension. The stochastic dynamics and the statistical properties of the 
polymer extension can be examined in a way very similar to the one developed in this 
letter and this will be the subject of separate publication. 

Our theory is based on the simple dumb-bell- like equation [jO.lJl . This equation is 
obviously approximate, taking into account only one variable (end-to-end vector R) 
of generally more complex dynamics. Therefore, it should be important to assess the 
effects of more realistic modeling. (This treatement should also account for internal 
conformational degrees of freedom.) 

Note that polymer tumbling was first observed in the steady shear flow experiments of 



Smith et al. (1999) and Hur et al. (2001) in which orientational fluctuations were driven 
by thermal noise, while our analysis has focused primarily on the case of tumbling driven 
by velocity fluctuations. Therefore, even though all of our results are directly appli- 
cable to the elastic turbulence setting of Groisman & Steinberg (2000,2001,2004), the 
examination of the statistics of the Langevin driven tumbling and angular distribution 
is, actually, a separate task. In particular, special attention should be given to the case 
of extremely elongated polymers in which the non-linearity of the polymer elasticity is 
essential. These Langevin-related problems will be examined elsewhere. Here, let us 
only note the universal nature of the exponential, large r, tail of the tumbling time PDF. 
Indeed, it is straightforward to check that the arguments presented above are generic, 
thus guaranteeing that, very much like in the velocity fluctuation driven case, the tail is 
also exponential in the Langevin-driven case. 
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